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We propose a mechanism for detection of a single spin center on a non-magnetic substrate. In the 
detection scheme, the STM tunnel current is correlated with the spin orientation. In the presence 
of magnetic field, the spin precesses and the tunnel current is modulated at the Larmor frequency. 
The mechanism relies on the effective spin-orbit interaction between the injected unpolarized STM 
current and the local spin center, which leads to the nodal structure of the spatial signal profile. 
Based on the proposed mechanism, the strongest spin-related signal can be expected for the systems 
with large spin-orbit coupling and low carrier concentration. 



There is no fundamental principle that precludes the 
single spin measurement. Possibility of a single spin 
observation is therefore a question of spatial and tem- 
poral resolution. The standard electron spin measure- 
ment technique - electron spin resonance - is limited to 
a macroscopic number of electron spins - 10 10 or more ; 
the state-of-the-art magnetic resonance force microscopy 
has recently achieved the resolution of about 100 fully 
polarized electron spins [^j. It has already been shown 
that optically induced ESR of a single spin is possible 
H . There are proposals for the spin detection using sin- 
gle electron transistor and spin-polarized current [§,1). 
In this letter we propose the theoretical basis for the new 
spin-detection technique - electron spin precession scan- 
ning tunneling microscopy (ESP-STM) - capable of sin- 
gle spin detection. 

The applications of single spin detection and manip- 
ulation range from the study of strongly correlated sys- 
tems, to nanotechnology, to quantum information pro- 
cessing. In the strongly correlated systems, the ability to 
detect single spin will allow one to investigate the mag- 
netism on the nano scale by detecting the changes in the 
spin behavior while entering magnetically ordered states 
H . One could also explore temperature evolution of the 
magnetic properties of a single paramagnetic atom in the 
Kondo regime [Q . Magnetic properties of spin centers in 
superconductors is another area where single spin plays 
an important role by generating intragap impurity states 
P,pJ . In nanotechnology, spins can be used as elementary 
information storage units. In the realm of quantum com- 
puting, several architecture proposals rely on the ability 
to manipulate and detect single spins |]To| , |Tl|] . Here we 
argue that STM offers a powerful technique to detect a 
single spin. 

The electron spin can be detected through its coupling 
to the magnetic field, which can either freeze the spin 
along the field, or make it precess around the field di- 
rection. To detect the single-spin signal, the superb spa- 
tial and temporal resolutions are required. The neces- 
sary level of sensitivity can be achieved in the Scanning 
Tunneling Microscope (STM). It is well known that spa- 
tial resolution of the STM is in the sub- Angstrom range. 
Therefore, the technique naturally lends itself for the de- 
tection of the single spin. In the STM setup, the spin 



precession in an external magnetic field can be detected 
through an ac modulation of the tunnel current. In this 
purely dc configuration, the ac tunnel current at the pre- 
cession, or Larmor, frequency can be generated due to 
the effective coupling of the precessing localized spin to 
the tunneling electrons. Hence we argue that STM can 
operate as Larmor frequency generator. The experimen- 
tal setup we consider is given in Fig. ^. 

Our work is motivated by experiments of Manassen et 
al. |Q. In these experiments, STM was used to mea- 
sure the tunneling current while scanning the surface of 
Si in the vicinity of a local spin impurity (Fe cluster) or 
imperfection (oxygen vacancy in Si-O). A small signal in 
the current power spectrum at the Larmor frequency was 
detected. The ac signal was spatially localized at the dis- 
tances on the order of 5 — 10 A from the spin site. The 
extreme localization of the signal and the linear scaling 
of its frequency with the magnetic field prompted Man- 
assen et al. to attribute the detected ac signal to the 
Larmor precession of a single spin site. 

Rather than going into further details of these experi- 
ments, we pose here a general question: Under what con- 
ditions it is possible to detect the single precessing spin 
with STM7 Surprisingly, we find that a number of mecha- 
nisms can generate a detectable ac current. The interac- 
tion that provides the coupling of the precessing spin to 
the tunneling current is the spin-orbit (SO) interaction. 
Based on our analysis, we argue that single spin detection 
is within the reach of current STM technology and can 
be realized in a variety of non-magnetic materials with 
isolated paramagnetic centers. 

Consider a localized magnetic site with spin S, S — 
1/2. In the presence of magnetic field, B, the spin-up 
(_E|) and spin-down (E^) energy levels are Zeeman split. 
At zero temperature, only the lowest energy state (par- 
allel to magnetic field) is occupied; at a finite tempera- 
ture, or due to an external excitation, the spin may be 
driven into the mixed state characterized by the wave- 
function, \tjj) = a\ j) + [3\ I), with the time evolution 
a = \a\ exp(— iE^t) and (3 — \f3\exp(—iEit + i(f)(t)). 
The drifting phase <p(t) determines the spin coherence 
time T<f, and is related to the ESR |T^| spin relaxation 
time Ti. In such a state, the spin expectation value, 
(4>(t)\S\tp(t)) / (il;(t)\ip(t)) , will precesses around the di- 
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rection of the magnetic field at the Larmor frequency, 



Hlot 



9^bB, 



(1) 



where g is the gyromagnetic ratio and /j,b is the Bohr 
magneton. In magnetic field of 100 Gauss, this frequency 
for a free electrons is 280 MHz. We assume hereafter that 
in the course of the ESP-STM experiment the spin is not 
in the groundstate either due to a thermal fluctuation or 
an external excitation by the tunneling current. We also 
assume that is sufficiently long for the spin preces- 
sion to be well defined, uj^t^ ^ 1. The validity of these 
assumptions is verified at the end of the paper. 
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FIG. 1. Experimental setup for the electron spin precession 
STM. In the applied magnetic field B, spin of the magnetic 
atom (e.g. Gd, shown in gold) is precessing around the field 
direction. When precisely positioned next to the spin site 
(within a few A), the STM tip can pick up an ac modulation 
of the tunnel current. 



We present now a mechanism for the ac STM current 
generation with no ac input based on the time-dependent 
modification of the tunneling density of states induced 
near the precessing spin in the presence of applied cur- 
rent. On the time scale of all conduction-electron pro- 
cesses, the spin precession is very slow. This follows 
from the comparison of the energy associated with the 
rf spin precession, Hlol ~ 10~ 6 eV, and the typical elec- 
tronic energy scale in metals or semiconductors, which is 
on the order of 1 eV. Hence, for every "instantaneous" 
spin orientation the electronic problem can be solved as 
if the local moment is static. This is similar to the Born- 
Oppenheimer approximation for atoms where atomic nu- 
clear dynamics is much slower than the electronic motion. 
We describe the system of electrons interacting with the 
local impurity spin by the Hamiltonian 



H = H a + JS • cr(0), 



(2) 



where J is the strength of the exchange interaction be- 
tween the local spin and conduction electron spin density, 
<t(0) = <t Q/3 c] 1 (0)c^(0), on the impurity site. To provide 
the coupling of the spin orientation to the orbital degrees 
of freedom of electrons, the non-interacting part of the 
Hamiltonian, Hq, should include spin-orbit interaction. 
We will focus here on the two-dimensional band of states 
on the material surface. From the symmetry considera- 
tions, the energy of the surface states contains spin-orbit 



part that is linear both in the electron spin, a, and mo- 
mentum, k, Jl4[ 



e Q /3(k) 



2m* 
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(3) 



Here, m* is the band mass of the surface states and gso 
is the SO coupling strength. The cross-product of the 
spin and momentum is projected onto the normal to the 
surface, n. 

The problem specified by equations @ and (||) can 
be solved for each instantaneous value of the precessing 
spin S(t). Specifically, we are looking for the correction 
to the conduction electron density of states due to the lo- 
cal spin. The Green function matrix for the conduction 
electrons is 



(4) 



In the presence of the current flow j in the system, 
the equilibrium momentum distribution is shifted by an 
amount proportional to the current, ko = jm*/ne, where 
n is the carrier density and e is electron charge. This 
shift can be approximately introduced by modifying the 
electronic Green functions fllal, 



G(k,w) 



(k-k ) 2 
uj gso k x er 

2m* ybUl 1 



(5) 



Here, only the kinetic energy momentum is shifted since 
SO term is explicitly non-invariant under shift (k — ► 
k k ) as it measures the SO effects due to lattice. The 
correction to the density of states near the Fermi surface 
caused by the impurity can be then calculated straight- 
forwardly in the first order in the scattering potential, 



6N(r) = ilm^G Q/3 (r)JS • a^G ia {-v) 



(6) 



ck/37 



Expanding in the strength of the SO interaction relative 
to the Fermi energy, the S-dependent contribution to the 
density of the surface states obtains in the first order in 
the exchange coupling and the strength of SO interaction, 



8N dN - 

■jf = gsoJ-^J ( k Fr)[k x SJ fi . 



(7) 



This correction depends on the distance from the spin 
center, r, through the Bessel function of the first kind, 
Jo (a;). This result also follows from considering Friedel 
spin oscillation of electrons coupled to the tunneling DOS 
in the presence of the SO term jl6| . It is time-dependent 
in the presence of magnetic field since the projection of S 
oscillates at the Larmor frequency jl7| . This expression is 
consistent with general properties of DOS. Since N(r, t) 
is a scalar, it should be invariant under time reversal, 
whereas S is odd under time reversal. Hence 5N(r,t) 
can depend only on product of the spin vector with some 
other vector that is odd under time reversal or on time 
derivative of spin vector. Possible combinations include 
SN ~ [I x S]n (Eq. |) and SN ~ [d t S(t) X r] ft . Indeed, we 
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have found mechanism for the latter contribution as well; 
however, for the experimental conditions of Refs. p2| the 
time-dependent DOS mechanism described here domi- 
nates. 

The current in the STM tunnel junction is proportional 
to the single electron density of states in the substrate. 
Therefore, the ac current contribution relative to the base 
dc STM current, /, can be estimated as 



SI{t) SN(t) 
I N~ ' 



(8) 



Notice, that the magnitude of the effect is proportional 
to the current in the system via ko. Experimentally, the 
necessary bias current can either be provided externally 
(by extra leads) or can be injected by STM itself. In the 
latter case, the direction of the current will vary in space 
and, in the absence of anisotropy, will flow directly away 
from the tip-sample junction. From the symmetry con- 
siderations, in Eq. (R), ko should be oriented from the tip 
towards the spin center. The magnitude of the equilib- 
rium shift is the 2D case considered here is fco ~ m*I/nr, 
where n is the 2D density of the surface carriers. Since ko 
itself is proportional to injected current, the S-dependent 
correction to the current is proportional to I 2 . For di- 
rect tunneling into the spin center, ko is parallel to the 
tunneling axis and hence there is no S-dependent contri- 
bution to the current. This nodal structure is a conse- 
quence of the problem symmetry. The magnitude of the 
time-dependent density of states is independent of the 
magnetic field orientation. 

The intuitive understanding of this result can be ob- 
tained using an analogue of the spin-dependent Biittiker- 
Landauer formalizm p8 19 1. The electrical circuit that 



corresponds to the STM experiment is shown in Fig. 
|2|. When current is injected through the tunnel junc- 
tion with resistance Rt, it can either flow to the leads 
directly (left branch with conductance T3) or through 
the channel affected by the spin center (right branch, 
and HQ. We assume here that there are no spin flips. 
From geometric considerations, the closer tip is laterally 
to the spin center, the larger is the ratio, r^)/!^. From 
the form of SO interaction, Eq. (^), the relevant spin a 
quantization axis for electrons flowing from the tunnel- 
ing point to the impurity is perpendicular both to the 
surface normal, n, and the vector connecting these two 
points, ko. The essence of the right branch of the cir- 
cuit is that in the presence of SO interaction the current 
becomes spin polarized, t\ ^ T\, and the transmission 
of the local magnetic center is spin sensitive, r| 7^ T^. 
Approximately, 



T° = T(l+g so k a/E F ), 
Y a 2 =T(1 + JSa/Ep), 



(9) 
(10) 



where Ep is the Fermi energy of the substrate, and S 
is the projection of spin S onto the a quantization axis 
(above). Similar to the Green function treatment pre- 
sented above, the expression for goes beyond linear 
response since it depends on the current through fco. The 



correction to the current due to the spin-selectivity of the 
circuit is 



<5J R 2 M Y gsokpJ 
I Rm + Rt Ep 



S, 



(11) 



where Rm = V(r + ^3). This expression is the same as 
the previously obtained result, Eqs. (0) and (g). 
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FIG. 2. Electrical circuit representing ESP-STM exper- 
iment. The current is injected through the tunnel junc- 
tion of resistance Rt- Inside the material the current 
can either flow through the impurity site (right, "impurity 
channel") or directly to the leads (left). In the impurity 
channel, S-dependence of the total current obtains through 
spin-current polarization due to spin-orbit coupling (conduc- 
tance rj Ij) and spin-dependent transmission through the 
impurity site (rj 7^ Tj). The resulting S-dependent contri- 
bution to the total circuit current, Eq. (0), is in agreement 
with the more sophisticated calculation, Eq. ((?]). 



From Eq. (0) it is clear that the largest ESP-STM 
effect can be obtained in the systems with large spin 
orbit interaction and strong particle-hole asymmetry in 
the density of states. Small density of surface carriers, 
n, helps to drive the carrier system out of equilibrium, 
which is crucial for the ESP-STM effect. All these pa- 
rameters can be tuned by the choice of the substrate ma- 
terial, as well as the type of the spin center, to achieve 
large response. For model parameters, gso^F ~ 0.01 eV, 
J ~ 1 eV, dN/dE ~ 1 cV~ 2 , and Fermi surface dis- 
placement ko/kp ~ 0.1, the effect in the density of states 
is 6N/N ~ 10~ 3 , which should be experimentally de- 
tectable. 

In the above discussion we assumed that the spin de- 
phasing time, T0, is long compared to the Larmor period, 
which implies weak coupling to environment. However, 
even under ideal circumstances, the measurement will in- 
duce spin dephasing. This backaction is an unavoidable 
side effect of measurement. We will now estimate the 
backaction of the tunneling electrons on the precessing 
spin. By applying a canonical transformation, the orig- 
inal Hamiltonian, Eqs. @ and (^), can be mapped to a 
model of two reservoirs (L and R) of spinless electrons 
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coupled through a spin-dependent tunneling matrix ele- 
ment @, 

H cS = H L + H R + h z S z + Yl 4 t T " + T i S xW- 

IGL, reR 

Here To is the tunneling amplitude that doesn't involve 
local spin and T\ ~ g so koJ(dN/dE)To is the spin- 
dependent part of the tunneling. The x axis is chosen 
along n x ko direction. For spin- 1/2 quantum systems, 
this model has been extensively studied [pl|j22j] . Two 
important parameters in the model are the electron tun- 
neling rate, l/r e = hT^NlNrV, and the measurement- 
induced spin-dephasing rate, l/r s ~ itTiNlNrV . It has 
been demonstrated p2| that under the weak measurement 
conditions, H/t s < h z , the stationary power spectrum of 
the tunnel current has a peak at the Larmor frequency, 



4/r s 2 



2e/ (w - oj L y + 1/t2 



For weak measurement, increasing applied voltage does 
not change the peak value relative to the shot noise level; 
however it increases the weight under the resonance peak. 
Under the stationary conditions, the coupling to the cur- 
rent fluctuations also provides the excitation for the spin. 
For the same parameters as above we find spin dephas- 
ing time l/r s ~ 10 4 . . . 10 5 Hz <C lo^. Experimentally, 
the observed linewidth is narrow, ~ 500 KHz <C lol [0 , 
which is consistent with our effective model. Any extrin- 
sic (unrelated to measurement) dephasing mechanisms 
will decrease the signal to shot noise ratio. 

The applicability of the presented ESP-STM mecha- 
nism is not limited to semiconductors. Other systems 
with stronger spin-orbit interaction, such as Au or Cu 
surfaces Q with local magnetic defects may be expected 
to show a similar effect. Yet another route is to use heavy 
magnetic atoms, such as Eu or Gd as the local spin im- 
purities with large intrinsic SO interactions. In this case, 
a result similar to Eq. |t] obtains, with the strength of 
the effect proportional to the spin-orbit coupling on the 
impurity site. In general, any form of SO interaction is 
capable of coupling the precessing local spin to the itin- 
erant electron density of states, which can be measured 
with STM. 

The spin-current coupling mechanism that we de- 
scribed applies to the static (dc) case as well. The dc case 
can be realized by applying a strong polarizing magnetic 
field and detecting the variation of the tunneling current 
as a function of the magnetic field orientation. This, 
however, is only possible at low temperatures, due to 
the limited strength of available magnetic fields. The dc 
measurement is further complicated by the low-frequency 
(1/f) noise. 

In conclusion, we have proposed a mechanism for the 
novel technique for the single spin detection with STM. 
It is based on the ac modulation of the tunneling cur- 
rent caused by the effective spin-orbit coupling to the 
precessing local spin. From the analysis of the coupling 



mechanism we identified classes of materials and experi- 
mental conditions for which significant single spin signal 
can be expected. 

We are grateful to G. Boebinger, F. Bronold, G. 
Brown, J.C. Davis, D. Eigler, C. Hammel, M. Hawley, 
Y. Manassen, A. Migliory, D. Mozyrsky, Y. Nazarov, 
D. Pines and A. Yazdani for useful discussions. We are 
grateful to D. Mozyrsky for pointing out unversal features 
of the current spectral density. This work was supported 
by the U.S. DoE. 



[1] M. Farle, Rep. Prog. Phys. 61, 755 (1998). 
[2] K. J. Bruland et al, Appl. Phys. Lett. 73, 3159 (1998). 
3 J. Kohler, J.A.J.M. Disselhorst, M.C.J.M. Donckers, 

(12) 

E.J.J. Groenen, J. Schmidt and W.E. Moerner, Nature 
363, 242 (1993); P. Tamarat, A. Maali, B. Lounis, and 
M. Orrit, J. Phys. Chem A 104, 1 (2000). 
[4] ] B.E. Kane, N.S. McAlpine, A.S. Dzurak, R.G. Clark, 
G. J. Milburn, He Bi Sun, and Howard Wiseman, Phys. 
Rev. bf B 61, 2961 (2000). 
[5] P. Recher, E. V. Sukhorukov, D. Loss, Phys. Rev. Lett. 

85, 1962 (2000). 
[6] S. Heinze, M. Bode, A. Kubezka, O. Pietzch, X. Nie, S. 
Blugel, and R. Wiesendanger, Science 288, 1805, (2000) 
and referencies therein; R. Wiesendanger, H.-J. Gun- 
terodt, G, Guntherod, R.J. Gamboni, and R. Ruf, Phys. 
Rev. Lett. 65, 247,(1990). 
[7] H. C. Manoharan, C. P. Lutz, and D. M. Eigler, Nature 

403, 512 (2000). 
[8] M.I. Salkola, A.V. Balatsky and J.R. Schrieffer, Phys. 

Rev. B 55, 12648 (1997) and references therein. 
[9] A. Yazdani, B. A. Jones, C. P. Lutz, M. F. Crommie, and 
D. M. Eigler, Science 275, 1767 (1997). 
[10] B. E. Kane, Nature 393 133 (1998). 
[11] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57 120 
(1998). 

[12] Y. Manassen, I. Mukhopadhyay and N. Ramesh Rao, 
Phys. Rev. B 61, 16223 (2000); Y. Manassen, R.J. 
Hamers, J.E. Demuth and A.J. Castellano Jr, Phys. Rev. 
Lett 62, 2531 (1989); D. Shachal and Y. Manssen, Phys. 
Rev. B 46, 4795 (1992); Y. Manassen, Journ. of Magnetic 
Reson. 126, 133 (1997). While working on this paper we 
became aware of a recent ESP-STM experiment by Cam- 
bridge group: C. Durkan et al., Appl. Phys. Lett., 80, 
458,(2002); H. Manoharan, Nature, 416, 24,(2002). 
[13] In semiconductors, the width of the ESR line I/T2 is 
quite narrow even at high temperatures, on the scale of 
10" 7 . . . 10" 8 s at T = 70 K in Fe-Si 2 ceramics. For phos- 
phoros in silicon, the spin relaxation time can be as long 
as 100 s. T. Miki, Y. Matsui, K. Matsubara and K. Kishi- 
moto, J. Appl. Phys. 75, 1693 (1994). See also R. Worner 
and O.F. Schirmer, Phys. Rev. B 34, 1381 (1986). 
[14] Yu. A. Bychkov and E. I. Rashba, JETP Lett. 39, 78 
(1984). 

[15] P. Lipavsky, V. Spicka, and B. Velicky, Phys. Rev. B 34, 

6933 (1986). 
[16] A.V. Balatsky, unpublished. 

[17] The effect of the magnetic field on electrons in the con- 



4 



LA-UR-01-XXXX 



February 6, 2008 



ventional metallic STM tip can be neglected due to the 
vanishingly small magnitude of the Zeeman energy rela- 
tive to the Fermi energy; the effect on the sample surface 
carriers can be neglected as well if the external magnetic 
field is relatively weak, such that [ibB < gsok.F- Fhis 
condition is satisfied in the systms with strong spin-orbit 
coupling that we are interested in. 



[18] M. Buttiker, Phys. Rev. Lett. 57, 1761 (1986). 
[19] H.X. Tang et ai, Phys. Rev. B 61, 4437 (2000). 
[20] I. Martin, unpublished. 

[21] S.A. Gurvitz, Phys. Rev. B 56, 15215 (1997). 
[22] A.N. Korotkov and D.V. Averin, Phys. Rev. B 64, 165310 
(2001) 



5 



